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ABSTRACT 

We make predictions for diffuse stellar mass fractions in dark matter halos from the scales of small 
spiral galaxies to those of large galaxy clusters. We use an extensively-tested analytic model for 
subhalo infall and evolution and empirical constraints from galaxy survey data to set the stellar mass 
in each accreted subhalo, which is added to the diffuse light as subhalos become disrupted due to 
interactions within their hosts. We predict that the stellar mass fraction in diffuse, intrahalo light 
should rise on average from ~ 0.5% to ~ 20% from small galaxy halos (~ 10 1:l Mq) to poor groups 
(~ 10 13 M Q ). The trend with mass flattens considerably beyond the group scale, increasing weakly 
from a fraction of ~ 20% in poor galaxy clusters (~ 1O 14 M ) to roughly ~ 30% in massive clusters 
(~ 10 15 Af Q ). The mass-dependent diffuse light fraction is governed primarily by the empirical fact 
that the mass-to- light ratio in galaxy halos must vary as a function of halo mass. Galaxy halos have 
little diffuse light because they accrete most of their mass in small subhalos that themselves have high 
mass-to-light ratios; stellar halos around galaxies are built primarily from disrupted dwarf-irregular- 
type galaxies with M* ~ 10 8 ' 5 M Q . The diffuse light in group and cluster halos is built from satellite 
galaxies that form stars efficiently; intracluster light is dominated by material liberated from massive 
galaxies with Af* ~ I0 11 M Q . Our results are consistent with existing observations spanning the 
galaxy, group, and cluster scale; however, they can be tested more rigorously in future deep surveys. 
Subject headings: Cosmology: theory - galaxies: formation - galaxies: evolution - clusters: diffuse 
light 



1. INTRODUCTION 

When Zwicky first observed the diffuse, luminous com- 
ponent of the Coma cluster of galaxies, i t was not clea r 
what processes were responsible for it (Z wickvl Il951| ). 
Today, the prevailing paradigm for structure forma- 
tion is hierarchical; galaxies and clusters of galaxies of 
all sizes are built through sequential mergers of many 
smaller objects. Hierarchical structure formation theo- 
ries provide a mechanism for the formation of intraclus- 
ter light as material lost fro m shredded galaxies over the 
course of clus t er formation (Gallagher & Ostrikc r Il972t 
Merrittj 119831 iBvrd fc Valtonenl ll99(TTD ubmski et all 
20031: iGnedirj 120031: iMihosI 1200 ' 



Murante et al. 2004 



Lin fc Mohr 2004; Wil lman et al . 2004; Sommer Larsenl 
20061 : iRudick et al.l 120061 : IConrov et alll2007f ). Whereas 
the building blocks of clusters are galaxies, galaxy- 
sized objects build their masses by acquiring rela- 
tively low-luminosity (dwarf) galaxies which may sub- 
sequently be destroyed by tides and heating processes 
to produce the dif fuse, stellar halos around galaxies 
like the Milky Way (fSearle fc Zirmlll978l; [Johnston et all 
1996HJohnstonllI998 : Bullock et al. 2 001 1 Johnston et all 
200ll: iBullock fc Johnston! 120051: [Robertson et al.l 1 2005 1 : 
Diemand et alj 120051 : iRead et alj|2006t iFont et alj|2006t 
Abadi et al.l l2006h ~Whether in clusters or galaxies, we 
refer to this diffuse material as "intra-halo light" (IHL) 
and adopt the symbol /ihl to express the fraction of the 
total system luminosity found in this diffuse component. 
In this paper, we explore the connection between the size 
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of a system and the relative fraction of its total light con- 
tributed by intrahalo stellar material. In particular, we 
predict the mean and variance in the IHL fraction as a 
function of dark matter halo mass, and we explore the 
origin of the scatter in IHL at fixed halo mass. 

Most of our knowledge about the IHL on galaxy 
scales (~ 10 11 — IO 12 M0) comes from star counts 
within the Local Group. The stellar halo of the 
Milky Way contains fiHL ~ 1% of the Galaxy's total 
luminosity (iMorrisonl 1993: iWetterer fc McGrawl H996I: 
Morrison et all 120001: IChiba fc Beers! 120001: lYa nnv et all 
2000tllvezic et alJl2000t ISiegel et al.ll2002t ). This number 
can be as large as /ihl ~ 2% if the unbound Sagittar- 
ius st ream stars are in cluded in the diffuse component 
{e.g., lLaw et al.ll2005l ). 

Interestingly, while the dark halo of M31 is thought 
to be roughly the sam e size as that of th e Milky Way 
(M M 3i ~ IO 12 M , see iKlypin et al.l l2002t ISeigar et all 
l2006a[ ). the recently-discovered, metal-poor stellar halo 
of M3I may contain a significantly higher fraction of 

that galaxy's light, /ihl ~ 2 5 — 5% (llrwin et al.l 

20051: iGuhathakurta et all 120051: iKalirai et alj 120061 : 
Chapman et al.ll2006l ). If the great Andromeda stream 
(jlbata et al.l l2001h were included as diffuse light, this 
count would be larger. These observations immediately 
suggest that there should be a substantial spread in IHL 
components among galaxy-sized systems. Detections 
of a stellar halo component in the smaller disk galaxy 
M33 (Mm33 ~ 10 n M w ) have recently b een reported 
(|Hood et alil2007tlMcConnachie et al.ll2006l ). These esti- 
mates are consistent with a very low stellar mass fraction 
in the M33 halo, /ihl 1%, although a higher number 
is not ruled out (A. Ferguson, private communication). 

In more distant galaxy halos, the IHL is both harder to 
detect and more difficult to discriminate from other ex- 
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tended components {e.g., iDalcanton fc Bernst ein 2002). 
Some results suggest that galactic stellar halos with 
/ihl ~ 1 — 5% are n o t uncommon (ISackett et al.l fl994l: 
Morrison et all Il997t IWeil et all Il997t iLequeux et all 
19981: lAbe et all ll999UZibetti fc Ferguson! \20oW. ReT 



cent work bv lBuehler et alJ (|2007t ) regarding the edge-on 
galaxy NGC 4244 indicates the existence of an asymmet- 
ric stellar component far above the system's exponential 
thin disk, although non-det ections are also reported in 
galaxies of similar size (e.g., IZheng et al.lll999tlFrv et all 
1999). Of particular interest is the case of NGC 300, a 
low-luminosity, late-type galaxy in which no stellar halo 
has yet been detected, despite the successful identifica- 
tion of an exponential disk that extends over 10 scale 
lengt hs from the disk's center (Bland-Hawthor n et al.l 
2005). Of course, the differences from object to object 
may reflect systematic observational issues, but taken at 
face value, they indicate that the IHL fraction around 
galaxy halos shows significant variation and that there 
may be a trend for low /ihl levels in small galaxies. 
Relevant determinations will become more precise as 
resolv ed-star surveys exte nd beyond the Local Group 
(e.g., Ide Jong et al1l2007| ). 

Diffuse light fractions on group scales (M v ; r ~ 
10 13 M Q ) also exhibit considerable variation from sys- 
tem to system; however, the IHL component typically 
accounts for a more substantial fraction of the total lu- 
minosity of the system than it does on galaxy scales. 
Observations suggest that the M81 and Leo groups have 
at most a few percent of their light in diff use form 
(|Feldmeier|[200l ICastro-Rodriguez et al.ll2003[) . At the 
opposite ext reme, HGC90 has a reported IHL fraction of 
/ihl ~ 45% (IWhite et al.l l2003). Studies in other groups 
of roughly the sa me size yield a range of IHL fractions , 
/thl ~ 5 - 30% (iDa Rocha fc Mendes de Oliveirall2005t 
lAguerri et al.l 120061 ). 

Galaxy clusters (~ 10 14 - 10 15 M Q ) typically show 
the highest fractions of diffuse light. Again the scat- 
ter in estimated values is s ignificant but values range 
from /ihl ~ 10 - 40% tiThuan fc Kormendvl \W7% 
Melnick et""aH Il977t lUson et alJ 119911: [Bernstein et al I 
19951: ICalcaneo-Roldan et al l l2000t iLin fc Mohrl 12004 



Feldmeier et alJ l2004at iMihos et al]|2005t iZibetti et all 



2001 iKrick et alJ [2 006: Seigar et al. 200gB). A review 
LJF 



bv lCiardullo et aLl (|2004h describes recent surveys in the 
small Fornax and Virgo clusters and points out a distinc- 
tive fall-off in the IHL fraction for systems smaller than 
L ~ 10 11 L Q - quite similar to the break we see in our 
predicted fractions below. Another interesting, though 
tentative trend, is that IHL fractions in clusters without 
cD galaxies appear to have a somewhat smaller typical 
/ihl(~ 10 — 20%), tha n do clusters with cD galaxies 
^Feldmeier et al.ll2004alrH ). 

Re cently, a series of papers by iGonzalez et all (|2007l 
l2005[) have argued that a more appropriate quantity to 
investigate is the sum of the diffuse intracluster light 
with that of the brightest cluster galaxy (more gener- 
ally, the "brightest halo galaxy" or BHG) since it is 
difficult to disentangle t he two components (the same 
appro ac h is advocated by IConrov. Wechsler. fc Kravtsovl 
2007). IGonzalez et all (|2005D find that the sum of 
IHL+BHG light is dominated by the diffuse component 
on cluster scale s IHL /(IHL + BHG) ~ 80%. Moreover, 
IGonzalez et alJ (|2007t ) find that, compared to the total 



light in the cluster, the IHL+BHG fraction decreases 
from ~ 35% in low-mass clusters M ~ 10 14 M Q , to ~ 25% 
in more massive clusters. As we discuss below, these 
trends are very much in line with our expectations. 

Comparing predictions for the IHL fraction with ob- 
servational data is a nontrivial task. On the galaxy 
scale, total stellar halo luminosities depend sensitively 
on the difficult-to-measure central core radius assigned 
to the faint halo component. In addition, the IHL will 
typically have a different color than the bound light in 
galaxies (because it likely traces different star formation 
epochs), implying that the IHL fraction should generally 
be a function of the luminosity band or tracer popula- 
tions used to determine it. Moreover, some traditional 
determinations of intracluster light have used relatively 
small patches of sky within the clusters themselves, in- 
troducing a statistical shot-noise error term into the in- 
ferred IHL values. The deep imaging necessary for in- 
tracluster observations is also heavily dependent on sky 
subtraction, providing another systematic barrier to pre- 
cision IHL measurements on these scales. Ideally, di- 
rect comparisons between predictions and observations 
will mimic the influence of particular observational tech- 
niques and choices on theoretical predictions. The goal 
of such studies would be to produce predictions and ob- 
servat ional results that can be compared in the ir detail 
(e.g., iRudick et all 120061 : lSommer-Larsenll2006fh 

In this paper, our aim is not to make such detailed com- 
parisons between predictions and observations. Rather, 
we focus on predicting the general behavior of IHL frac- 
tions as a function of the size of the system from dwarf 
galaxies to large clusters 4 . We also explore the typ- 
ical galaxy size that contributes to IHL as a function 
of halo mass and explore the scatter from system to 
system at fixed host mass. The scope of this study 
represents a challenge for direct numerical simulation 
of halo formation due to the limited dynamic range of 
such computations. To achieve our goa ls, we rely on 
an an alytic treatment of halo formation (|Zentner et al.l 
120051 see below). We normalize the stellar content of 
our accreting halos t o match empirical constraints from 

z 

l200l 

make the explicit assumption that stellar material in 
galaxies is liberated when their dark matter halos be- 
come significantly stripped. We make no distinction be- 
tween material that has recently been liberated by tidal 
interactions (which may therefore appear as stream-like 
structure) and the general diffuse background. In order 
to avoid any ambiguities associated with the evolution of 
luminosity in different components, we quote the diffuse 
stellar mass fraction, f mL = Mf ff /M* otal . 

In the next section, we outline our two-step model for 
IHL predictions. In § [31 we briefly describe a toy model 
for the scaling of the IHL fraction with halo mass that 
serves both to frame our expectations for the fiducial 
result and to demonstrate the generality of this scaling. 
In § 4, we present our results for IHL fractions, reserving 
§ [5] for discussion and review. Throughout this work we 



iccreting halos t o match empirical constraints from 
observations (lYang et alJ 120031: IVale fc Ostrikerl 
[ iBell fc de Jorigl 120011 : Ide Jong fc Belli l2006h . We 



4 Note that Murante ct al. (2004) predicted a positive trend be- 
tween intracluster light fraction and cluster mass, focusing only 
on cluster scales. We extend the range of mass by more than two 
orders of magnitude. 
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Fig. 1. — The differential mass fraction in subhalos of mass M sa t, d//d(logM sa t), as a function of satellite mass. We plot the mass 
fractions for host halos of five masses as marked in each panel. The solid lines represent the mass fraction of all satellites accreted throughout 
the entire history of the host system. The red dashed lines represent the fraction contributed by subhalos that are eventually disrupted 
according to our algorithm, while the green dot-dashed lines represent halos that survive according to our algorithm. The upper horizontal 
axes show the luminosities of the galaxies assigned to each subhalo. 



adopt a ACDM cosmology model with h — 0.7, fl m = 
1 — Q\ = 0.3, and a primordial power spectrum which is 
scale-invariant, n = 1, and normalized to erg = 0.9. 

2. METHODS 

2.1. Dark halo accretion and disruption 

We model host dark matter halo mass accretion his- 
tories and track the evolution of accreted dark matter 
subhalos using an analytic prescription developed and 
tested against dissipa tionless cosmological simulations by 
IZentner et al.l (|2005l Z05 he reafter). This approach i s 
based on the earlier model of IZentner fc Bullockl (|2003f ). 
The analytic technique enables us to explore quickly the 
expected variety of accretion and disruption histories for 
host halos at a series of different masses. The model 
has proven remarkably successful at reproducing subhalo 
count statistics, radial distributions, and two-point clus- 
tering statistics measured in full, high-resolution iV-body 
simulations in regimes where the two techniques are com- 
mensurable. This success spans more than 3 orders of 
magnitude in host halo mass and persists as a function 
of redshift (Z05). The range over which this agreement 



is known to exist is limited only by the dynamic range of 
the simulations used by Z05. In what follows, we apply 
the analytic model outside the range over which it is well 
tested, but we know of no reason that it should fail out- 
side of this range. Of course, more precise estimates that 
involve full A^-body and hydrodynamical simulations will 
need to be made to refine our predictions; however, the 
general success of the model suggests that our predictions 
should be accurate enough that the approximate dynam- 
ical treatment of subhalos is not the limiting source of 
error and that potential differences are likely to test our 
assumptions about the evolution of stellar mass. Even so, 
many of the qualitative trends we derive are reflections of 
very general features of hierarchical structure formation 
and should be robust. Here we provide a brief overview 
of the technique and refer the reader to IZentner et al.l 
(120051) and the similar models ofrTavlor fc Babull (12004) 
iPenarrubia k. Bensonl (120051). iFaltenbacher fc Mathewsl 
(|2005h . and Ivan den Bosch et alJ ( 20051 ) for more detail. 

In hierarchical cosmologies like ACDM, dark halos ac- 
cumulate their mass through a series of mergers with 
smaller objects. The first step in our model is to se- 
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lect a host halo mass Mh os t at z — and gener- 
ate a subhalo-based mass accretion history using the 



extended Press- Sch echter formalism (iBond et al 



lLacev &: Coldll993l for a recent review see lZentne" 20061 ) 



1991; 



We use the particu l ar im plementation advocated by 
iSomerville fc Kolattl (fl999h . The merger tree contains 
a list of all of the merger times and masses of all of the 
smaller halos that have merged to form the final object. 
Every time there is a merger, the smaller object becomes 
a subhalo of the larger object. This is a Monte Carlo pro- 
cedure. Each merger event is drawn from a probability 
distribution and by realizing merger trees for numerous 
halos of the same final mass, we can probe the variety of 
formation histories that lead to final objects of the same 
size. As we discuss below, this variety of halo mass ac- 
quisition histories is a primary source of scatter in the 
fraction of IHL at fixed host mass. 

After constructing a large number of merger histo- 
ries at each final mass scale, we then track the evo- 
lution of subhalos in the dense environments of their 
host systems. Specifically, we assign an initial orbital 
energy and impact parameter to each merging subhalo. 
These values are chosen from probability distributions 
extracted from cosmological A^-body simulations in Z05. 
We then integrate the orbit of each subhalo in the po- 
tential of the main halo from the time of accretion to 
the epoch of observation. We model tidal mass loss us- 
ing a modified tidal approximation and a prescription for 
internal heating, as well as the effect of dynamical fric- 
tion using an adaptati on of the Chand rasekhar formula 
Chandrasekharl I1943D suggested by lHashimoto et al.l 
2003). For simplicity, we model the density structures of 



all halos and subhalos by th e sphe rically-symmetric den- 
sity profile of iNavarro et al.l (|1997l NFW) . For each halo 
and each subhalo, we set the concentrati on of the NFW 
profile according to the prescription of Wcchslc r et al.l 
(2002) to account for the correlation between mass accre- 
tion history and halo concentration. Masses are defined 
relative to the virial overdensity A V j r , where A V i r = 337 
at z = (e.g., iBullock et aLll2001[ ). 

Each subhalo has a well-defined rotation curve V c = 
y/GM(< r)/r, that peaks at a velocity V ma , x . As the 
subhalo orbits within its host, it gradually loses mass at 
all radii and the value of V max declines. A subhalo is 
declared to be "disrupted" when its maximum circular 
velocity falls below V^ rit = f„itV max (tacc)- The quantity 
/ cr it is a parameter that allows us to determine when 
the galaxy associated with each halo will contribute its 
stars to the diffuse light of its host halo. We have some 
freedom to tune / cr ;t to match empirical constraints on 
the n umber of survivin g satellite galaxies per halo (see, 
e.g., lYang et al.ll2003l and discussions below). We ex- 
pect that a satellite galaxy will typically remain bound 
within its subhalo until the subhalo loses a significant 
portion of its mass. Physically, f cr i t should not be so 
high that a system would be classified as disrupted when 
its host halo is only slightly less massive than it was at 
accretion. Similarly, a very low choice of f cr it would en- 
sure that the galaxy would not be considered destroyed 
until the dark matter in its host subhalo is less massive 
than the galaxy itself. 

Adopting a simple mass-scaling argument may allow 
us to gain physical insight into the disruption threshold, 



if we consider that the v irial mass of a halo sc ales approx- 
imately as M oc V%£ x (jBullock et al.ll2001[ ). With this 
in mind, an f cr it value of 0.8 translates to the halo being 
"disrupted" when it has lost just over half its mass, while 
fcrit — 0.2 implies a mass-loss threshold of more than 
99.5%. Clearly, the smaller our f crit is, the more assured 
we can be that galaxies meeting the criterion are truly 
dispersed, but if this parameter is chosen to be too small 
then we may falsely associate galaxies with what should 
rightly be diffuse, luminous material. As discussed be- 
low, we adopt fcrit = 0.6 as our fiducial value primarily 
because it produces reasonable agreement with empiri- 
cal constraints described in § 12.21 This choice implies 
disruption begins to occur when just under ~ 20% of the 
halo mass remains bound. 

Our definition of "disruption" is not necessarily meant 
to indicate that beyond this threshold, a subhalo must 
become physically unbound due to the interaction within 
the host potential. Rather, our intention is to introduce 
some effective criteria whereby it would be sensible to 
assign a large fraction of the subhalo's stellar mass to 
a diffuse component. The parameter f cr n denotes this 
transition from a bound galaxy component that con- 
tributes little diffuse light, to a tenuous structure that 
relinquishes most of its stellar mass to the diffuse com- 
ponent of the host halo. In our IHL predictions, we make 
the explicit assumption that the stars initially assigned 
to a subhalo become "diffuse" when that subhalo is "dis- 
rupted" according to the aforementioned criterion. 

Armed with a prescription for the mass accretion histo- 
ries of halos and the subsequent orbital dynamics of their 
satellites, we can investigate the predicted substructure 
distributions and overall accretion histories for host halos 
of various masses. Our main results rely on 1000 real- 
izations for virial host masses from 10 10 5 to 1O 15 O M , 
with four discrete intervals in each decade (e.g., in log- 
space: 11.2, 11.5, 11.8, 12.0, etc.) for a total of 19 
mass bins. The solid lines in Figure [1] show the frac- 
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Fig. 2. — The total mass-to-central-galaxy-light ratio as a 
function of halo m ass. The solid curve is the value inferred by 
lYang et al.l fl2003fl . This represents the L C (M) relation that we 
adopt in our fiducial models. Fo r comparison, the dotte d line is the 
mass-to-light ratio presented by Vale & Ostriker (2004). 
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tion of host halo mass accreted in satellite halos of a 
given mass (d//dlog M sa t) averaged over 1000 realiza- 
tions for host halos of mass Mh os t = \0 11 Mq (top left) 
through Mhost = 10 15 M Q (bottom left) at z = 0. To 
be explicit, /(> M sat ) is the cumulative mass fraction 
in satellites larger than M sat and our prescription de- 
mands that /(> M sat ) — ► 1 as M sa t — > (i.e., all of a 
halo's mass is accreted in subhalos of some size). In each 
panel, the dot-dashed lines include only subhalos that 
survive to the present day and the dashed lines include 
only subhalos that are disrupted, according to the above 
definition, between the epoch of accretion and z = 0. 

It is important to note that regardless of host mass, 
the majority of mass is accreted in subhalos of mass 
M sa t ~ 0.05 — 0.1 Mhost- In addition, surviving subhalos 
contribute much less mass than their destroyed coun- 
terparts of similar size in galactic systems, while their 
relative contributions are more even in cluster-size halos. 
This trend arises because high-mass halos accrete their 
subhalos more recently than low-mass halos. Therefore, 
the subhalos of low-mass halos are typically more dy- 
namically evolved and more likely to be destroyed (see 
Z05). These facts are fundamental to understanding the 
diffuse light fractions as a function of host mass, the con- 
sequences of which we explore in § 3. 

2.2. Assigning Light to Dark Matter Halos 

We assign a luminous component to each accreted halo 
using an empirical model that is normalized to z = 
galaxy constraints. We assume that every accreted sub- 
halo and every host halo contains a central galaxy. For 
every system accreted at time i acc we determine the stel- 
lar mass that this system would have today (at t = to 
or z — 0) according to empirical mass-to-light ratios. 
Next, we extrapolate this z = value backward in time 
to obtain M*(f = t acc ) using an empirically-motivated 
star formation law. The z = normalization guaran- 
tees that our model produces the required relationship 
between host halo mass and (central) galaxy luminosity 
required to match local galaxy counts and galaxy clus- 
tering observations. 

As we show below, our results for IHL fractions are 
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Fig. 3. — The cumulative number of surviving galaxies as a 
function of luminosity, both in our fiducial model (solid) and in the 
Yang et al. CLF analysis (dashed), for a host halo of mass M^ oat = 
10 14 - 5 Mq. The model line represents the mean of 1000 fiducial 
realizations, and the error bars representing the error on the mean 
over the sample are smaller than the model line's thickness on this 
plot. 



quite insensitive to star formation assumptions. In- 
deed, our primary prediction, that the IHL fraction in 
halos will vary strongly with mass scale, is extremely 
robust, and is driven by the empirical fact that the 
global mass-to-light ratio (M/L) in ACDM halos must 
vary strongly with host halo mass in order to reproduce 
the observed gal axy luminosity funct i on and clustering 
statis t ics {e.g., IWhite fc Reesl 119781: iKauffmann et al.l 
19931 ISomerville fc Primackl Il999t iTinker et all 120051: 
Coorav fc Milosavlievicll2005l ). 

We adopt the M/L relat ion inferred by 
lYang. Mo. fc v aii don Bosch (2003) in their model 
"Ml." lYang et all ([20031 ) used data from the Two- 
Degree Field Galaxy Redshift Survey (2dFGRS) to 
constrain the "conditional luminosity function" (CLF) 
of the 2dFGRS galaxies. This comparison allowed them 
to derive a characteristic B-band luminosity, L C (M), for 
the central (brightest) galaxies tha t sit in halos of virial 
mass M (for related analyses, see [v an den Bosc h et "all 
2001 ITinker et alj 120051 fcoorav fc Milosavlievicl [20051 : 
Yang et all |2005| ) . The solid line in Figure [2] shows the 



inferred total mass-to-light ratio (M/L c ) as a function 
of halo mass 5 . The dotted li ne sh ows an independent 
result from IVale fc Ostrikerl (200 \ ). which we utilize 
below in order to investigate the dependence of our 
conclusions on the specific choice of (M/L c ) function. 
Note that in both cases, galaxy formation is most 
efficient in dark halos of virial mass M ~ 5 x 10 M© 
and the conversion of baryons to stars is increasingly 
less efficient as we consider halos with masses either 
larger or smaller than this scale. 

In practice, we work with stellar mass rather than 
luminosity to avoid uncertainties associated with stel- 
lar population evolution. After computing the central 
galaxy luminosity using the Yang et al. relation shown 
in Figure [51 we convert this luminosity to a stellar mass 
using the average "m ass-dependent dust" relation from 
iBell fc de Jonel pOOlh : 



Ml 
L c 



0.75 



0.33 



(1) 



We have adjusted the IBell fc de Jong] (|200lD normaliza- 
tion down by a factor of 1.26 as ad vocated by their more 
recent work ()de Jong fc Bell 2006]). In the final analysis, 
our predictions for diffuse light fractions depend very lit- 
tle on the overall normalization. 

If we were interested only in contemporary galaxy and 
halo properties, the (M/L c ) relation at z = would suf- 
fice. However, the majority of the subhalos in our mod- 
els are accreted well before z = 0. This fact forces us 
to adopt a star formation prescription in order to ex- 
trapolate our z — stellar masses to earlier times. For 
simplicity, we assume that a galaxy's star formation is 
truncated at the time it is accreted into a larger host, 
perhaps due to ram pressure stripping or the fact that 
gas leaks more readily out of the potential well of a sub- 
halo located in a background host than it would if the 
satellite were left alone in the intergalactic field. 

After setting the z = stellar mass, we adopt a sim- 
ple approach that models star formation with minimal 
parameterization, in order to estimate the stellar mass 

5 Data table kindly provided by X. Yang. 
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that a particular system would have had at the time of 
accretion, t acc < to ~ 13.6 Gyr. We impose a history 

'to -t' 



M,(t) = M*(i ) 



1 



to 



(2) 



This equation introduces a second free parameter a into 
our analysis, which can be adjusted to produce a wide 
range of evolutions for the stellar mass in a system. For 
example, a = 0.25 will cause a galaxy to form most of 
its stellar component within the last two Gyr, while a 
larger value of a — 2 results in a system with a much 
earlier formation epoch, increasing the lookback time to 
half-stellar-mass formation by roughly a factor of five. 
As in our choice of f cr it = 0-6 for the disruption pa- 
rameter, we similarly adopt a = 1 to best match the 
expected luminosity functi on of satellite gala xies in host 
halos of a given mass from lYang et a l. (2003). We make 
these choices primarily for convenience and concreteness, 
and we demonstrate in § 4.3 that our main results for 
IHL fractions are largely insensitive to these parameter 
choices. 

An example of our (surviving) galaxy population is de- 
scribed by the cumulative luminosity function plotted in 
Figure [3l We caution that this figure, unlike our main 
results below, focuses on galaxy luminosity rather than 
stellar mass. While we allow for stellar mass buildup 
with time, we do not include any luminosity evolution, 
which should be important for determining the B-band 
luminosity of cluster galaxies. We would expect, for ex- 
ample, that systems that have survived in the cluster 
environment for several Gyr would have stopped form- 
ing stars and faded in blue light. Instead, we have used 
Equation 1 to convert between stellar mass and lumi- 
nosity regardless of the redshift at which the satellite 
was accreted. We neglect any explicit stellar population 
modeling in order to keep our methods as simple as pos- 
sible and to concentrate on robust, model-independent 
predictions. We present this only to demonstrate the 
gross consistency with inferred satellite galaxy popula- 
tions in halos and do not adopt this strategy for any of 
our predictions below. 

The solid line in Figure[3]shows the cumulative number 
of surviving galaxies (including the central galaxy) in a 
cluster-sized host halo, Afhost = 10 14 ' 5 M Q , as a function 
of galaxy luminosity. The dashed line shows the empiri- 
cally derived CLF result of lYang et alj (|2003t ). Here, and 
for the rest of the paper unless otherwise stated, we have 
used our fiducial parameter choices f cr u — 0.6 and a = 1. 
Overall, the agreement is encouraging, and we find simi- 
lar results for host halos of various masses. We match the 
empirical expectation quite well for the brightest galax- 
ies, which is not surprising because the central galaxy is 
forced to be of the "correct" luminosity by construction. 
We gradually begin to over-predict satellite galaxy counts 
relative to the empirical line at faint luminosities, but as 
we now argue, this is not of serious concern for a number 
of reasons. First, as we show below, the vast majority 
of accreted stellar mass will be contributed by the most 
massive accreted galaxies. This suggests that an accu- 
rate reproduction of the brightest satellites is the most 
important aspect of the IHL calculation. Second, the 
faintest galaxies will likely be most affected by luminos- 
ity evolution (which we do not include). These objects 
tend to survive the tug of dynamical friction longer than 



their more massive companions, and we expect them to 
fade considerably in B-band light as they evolve in the 
cluster environment. Finally, though errors in the de- 
rived luminosity f unction are not explicitly discussed in 
I Yang et"aTI (|2003h . the faintest galaxies in clusters are 
certainly weakly constrained by gross galaxy statistics 
because they are only a minor contributor to the global 
count of faint galaxies in t he universe (see, e .g., the clus- 
ter luminosity functions in I Yang et al.|[2005T ). 

2.3. Evolving the Diffuse Stellar Mass and the Central 
Galaxy Stellar Mass 

To calculate the amount of diffuse light in a clus- 
ter, group, or galaxy halo, it is necessary to determine 
whether the stellar material from a disrupted halo should 
be included as extended, diffuse material or as material 
that is incorporated into the central galaxy. In practice, 
infalling satellites should deposit stellar mass into both 
the central galaxy and the diffuse component. However, 
modeling these interactions in detail is challenging, so 
it is difficult to budget the fraction of the infalling stel- 
lar material that should be assigned to the diffuse com- 
ponent and the fraction that should be assigned to the 
central galaxy. 

To circumvent this complication, we employ two sim- 
ple, alternative models for adding stellar mass to the cen- 
tral galaxy and diffuse components that should bracket 
the outcome from a full modeling of the baryonic com- 
ponents. In Case 1, we classify all stellar material from 
disrupted subhalos as IHL. In this case, the diffuse stel- 
lar mass fractions should be maximized. In Case 2, we 
exclude from the IHL all galaxy stars from subhalos that 
make an approach closer than a radius r c (M^ ost ) to the 
center of their host halos. In these instances, we add 
the liberated stars to the stellar mass of the central ob- 
ject. Relative to Case 1, stellar mass is removed from the 
diffuse component and added to the light of the central 
galaxy. This causes diffuse stellar mass fractions to be 
smaller in this case. We associate r c with a characteris- 
tic outer radius for the central galaxy. To be conserva- 
tive, we adopt a fairly large outer radius r c — 10 kpc, 
for central galaxies of stellar mass M* = 4 x 10 10 M Q . 
We assume that r c scales according to the findings of 
IShen et all (|2003l ) for Petrosian half-light radii of galax- 
ies in the Sloan Digital Sky Survey. To be explicit, we 
use r c oc M, 0,4 for M* > M* and r c oc M" 16 otherwise. 
The generous value of r c along with the assumption that 
all stellar mass is assigned to the central galaxy and that 
none of the stellar mass goes into the diffuse component 
should lead to minimum IHL fractions in Case 2. 

Thus far, we have only considered the disruption of 
subhalos belonging to the trunk level of the host halo's 
merger tree, i.e., we have not made any determinations 
about the diffuse stellar content already present in ac- 
creting subhalos, often r eferred to as "pre-p rocessed" in- 
trahalo light (see, e.g., iRudick et al.1 120061 ). We do not 
expect galactic-scale host halos to carry much of this pre- 
processed stellar material, since accreting dwarf satel- 
lite galaxies typically have very little of their luminosity 
in diffuse form, but cluster-sized hosts accrete most of 
their mass in galaxy groups that may have a significant 
amount of IHL already present. In order to replicate this 
phenomenon, we first obtain our fiducial result (without 
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Fig. 4. — The diffuse light fraction as a function of host halo mass, for systems with virial mass between 10 la5 M Q and 10 15 M Q . In both 
panels the diamonds denote the mean of the distribution of IHL fractions at fixed mass based on 1000 realizations of our analytic model. 
The light shaded region shows the 95% range of the distribution of IHL fractions at fixed mass and the dark shaded region contains 68% of 
the distribution. The solid lines show the median of the distribution. Note that the median differs markedly from the mean at small host 
masses, illustrating the skewness of the IHL distribution in that range. The dashed line represents the preparatory IHL fraction, without 
the addition of pre-processed diffuse material already in subhalos at the time of accretion. The left panel shows results from Case 1 in 
which all stellar material from disrupted subhalos is included as diffuse light. The right panel shows Case 2 in which the subhalos on orbits 
passing within a critical radius r c of the host halo center contribute their light to the central galaxy rather than t he diffuse compone nt (see 
§[2). The upper axes show the corresponding central galaxy (BHG) luminosity derived from the prescription in I Yang et al.| (120031 ), with 
the Case 2 axis accounting for the average amount of stellar mass contributed by subhalos that merge with the BHG. 



the presence of pre-processed IHL), which is then used 
to interpolate an initial IHL value for each accreting 
subhalo. We then reproduce the fiducial IHL fraction, 
this time including diffuse stellar mass already present 
in subhalos and contributing that amount to the host's 
total IHL upon the subhalo's accretion, essentially boot- 
strapping Case 1 into itself in order to account for pre- 
processed intrahalo luminosity. We expect this model to 
differentiate itself from the initial result on large mass 
scales, at which subhalos are likely to have diffuse stel- 
lar components that contribute a non-negligible portion 
of the total subhalo luminosity. It is also worth noting 
that this second-order substructure is likely to carry in- 
creasingly less IHL than the fiducial model predicts at 
fixed mass because these subhalos will be younger and 
less luminous and will be less dynamically evolved so less 
luminosity could have been converted to diffuse luminos- 
ity. This indicates that our pre-processing method will 
slightly overestimate the contribution to the total host 
IHL made by diffuse stellar mass belonging to higher- 
order subhalos. 

3. A TOY MODEL FOR THE INTRAHALO LIGHT 
FRACTION 

Before proceeding, we derive a crude, analytic estimate 
for the scaling of the IHL fraction as a function of host 
halo mass. This model serves to frame our expectations 
for the general behavior of IHL fraction with mass, to 
highlight the features of hierarchical structure formation 
models most relevant to the determination of IHL frac- 
tions as a function of halo mass, and to demonstrate the 
generality of the halo mass-IHL fraction trends that we 
present in more detail in the following section. 

The gross scaling of IHL fraction with halo mass can be 
understood from two robust, cosmologically-motivated 
inputs: 

1. Host halos of mass M tend to accrete most of their 
mass in subhalos of mass M sat ~ 0.05 — 0.1 Afhost 



(Figure Q]) , and these halos are disrupted very effi- 
ciently due to dynamical friction. 

2. Galaxy formation picks out a typical halo mass 
M t ~ 5 x 10 n -Mo, where star formation is most 
efficient, and the efficiency of star formation de- 
clines rapidly away from this value (Figure [2|). 

To begin with, it is useful to introduce an approximate 
analytic fit to the adopted (M/L c ) relation from Yang et 
al. (2003): 
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The differential contribution to the IHL fraction from a 
satellite of mass M sat can be computed by introducing 
two parameters: /destroy that encapsulates the probabil- 
ity that this satellite will be destroyed, and /dig which 
describes the fraction of the satellite's stellar mass that 
contributes to the diffuse light once it is destroyed. Con- 
ceptually, this decomposition is useful, because /destroy 
has a known dependence upon host and satellite halo 
masses (Z05, see Fig. Q]). We will show that this de- 
pendence is subdominant, so for our purposes we can 
condense these into a single parameter /d = /destroy /diff 
that accounts for the average fraction of its total stellar 
mass that a satellite contributes to the IHL. As we stated 
above, the mass dependence of /destroy is weak and is not 
the dominant factor that gives rise to the mass scaling of 
the IHL fraction and, for simplicity, we will assume the 
composite parameter /d to be a slowly-varying function 
of mass. 

The differential contribution to the IHL fraction from 
satellite halos in the mass range dM sat around M sat is 
then 



d/iHL j. L(M sat ) dn acc 
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Fig. 5. — Similar to Figure[4] where now we plot the IHL fraction relative to the sum of the IHL and the brightest halo galaxy (BHG) . 
The rise at high host masses compared to Figure [4] reflects the fact that much of the stellar mass in clusters is bo und to satellite galax ies. 
Note that the IHL mass dominates that of the central galaxy mass on cluster scales, in accord with observations (Gonzalez et al. 2005). 
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where dn acc /dM sat is the mass function of accreted satel- 
lites. In general the amount of total and stellar mass 
accreted into the system is dominated by the few most 
massive satellites near ~ M host /20 (see Fig. [T]). As a 
final rough approximation, we assume that satellites of 
this mass dominate the integral over M sa t. This gives 



/iHL(M)~/ d n e ff 



L(M host /20) 
L(M } 



host ) 



-0.005/d n eS 



(M host /M t f 2 



1 



(M host /20M t ) 3/2 



(6) 



(7) 



where we have introduced a final parameter n er r which 
represents an effective number of satellites near mass 
Msat = A/host /20 and will be of order unity (Fig. [1} and 
/d is understood to be evaluated near M sat = Mhost/20. 

As will be clear in the following section, this extremely 
simple model captures the general features of our more 
detailed predictions. In our full model, /d should be less 
than one and n c g should be of order unity. This sim- 
ple model predicts that the IHL fraction should have a 
small and nearly constant value below M t ~ 5 x 10 11 M Q , 
/ihl ~ 5 x 10~ 3 . We expect a rapid rise in the IHL frac- 
tion with halo mass, /ihl oc M 3 ' 2 , for halos in the mass 
range M t < Mi lost < 20M t . In physical units this range is 
5xlO n M Q < M host < 10 13 M Q , and represents the range 
of transition between MW-like galaxies and small groups 
of galaxies. For host halos more massive than groups, 
Mhost >, 20M t ~ 1O 13 M , both the relevant satellite ha- 
los and host halos fall along the power-law regime of the 
Mj L c function and we expect the IHL fraction to remain 
roughly constant, /ihl IS 40%. 

At this point, it behooves us to summarize the points 
that this model illuminates regarding the IHL on differ- 
ent scales. In our model, it is approximately true that 
only the relative sizes of host and satellite objects deter- 
mine the probability for satellites to deposit their stellar 
mass into the diffuse component. Halos acquire most of 
their mass, dark or stellar, in a relatively small number 



of accreting objects of order 1/20 the size of the par- 
ent object (see Fig. Q] and Fig. [7] in the following sec- 
tion). Though the details are not known, it is an em- 
pirical fact that in a hierarchical cold dark matter cos- 
mology, the process of galaxy formation must pick out a 
halo mass scale where galaxy formation is most efficient 
(Al t ~ 5 x 10 11 M Q ) and that this efficiency drops at 
both lower and higher masses. Halos less massive than 
~ 20Af t will accrete little stellar mass in satellite objects 
and thus have little opportunity to build a diffuse, stellar 
halo. Halos more massive than ~ 20M t will accrete many 
satellite halos with masses such that they form stars at 
near peak efficiency. As these host halos bring in satel- 
lites with lots of stars, they have ample opportunity to 
build diffuse stellar halos. The general conclusion that 
diffuse light fractions should increase from very small val- 
ues in galaxy-sized systems to larger values in group- to 
cluster-sized systems seems difficult to avoid in the con- 
text of hierarchical cold dark matter structure formation. 

4. RESULTS 
4.1. IHL Fraction and Dark Halo Mass 

The two panels of Figure |4] show our primary re- 
sults. The predicted diffuse stellar mass fraction, /ihl = 
M^ lff /M' otal , is shown as a function of host halo mass. 
The total stellar mass M* otal includes the stellar mass in 
the diffuse component (IHL), satellite galaxies, and the 
central galaxy. In this section we will refer to the central 
galaxy as the "brightest halo galaxy" or BHG, in analogy 
with the brightest cluster galaxy (BCG) in clusters. The 
left panel shows results for Case 1 in which we assign all 
light from disrupted subhalos to the diffuse component. 
The right panel shows Case 2 in which we exclude any 
stars that were in subhalos having passed within r c of the 
host halo center from the diffuse component and instead 
add this material to the stellar mass of the central BHG. 
Diamonds show the average value of /ihl and the thin 
solid line shows the median. These results are derived 
from 1000 realizations for each host halo mass. The light 
and dark shaded regions span the 95% and 68% regions 
of the distribution respectively, centered on the median. 

The upper axis in the Case 1 (left) panel of Figure 0] 
shows the luminosity of the central galaxy according to 
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the I Yang et ail ((2003) mapping, while the Case 2 (right) 
panel upper axis label gives the mean BHG luminosity at 
a given host mass that we obtained by averaging the to- 
tal merged subhalo luminosity (plus the assigned central 
galaxy luminosity). We should note here that our Case 2 
model is not self-consistent, in that we first assign a cen- 
tral galaxy luminosity according to the z = conditional 
luminosity function and then subsequently add stellar 
mass via subhalo mergers, which will obviously produce 
incorrect present-day stellar mass functions. However, 
on cluster scales, the Case 2 central galaxy luminosi- 
ties are many tim es larger than that required by the 
lYang et alj (|2003f) analysis (e.g., L B hg ~ 9 x lO n L 
compared to ~ 2 x lO 11 ^© from Yang et al.), implying 
that even if our BHG were composed entirely of merged 
material (without any stellar mass produced by direct 
cooling processes), we would still over-predict the lumi- 
nosity of the central object. We therefore present the 
Case 2 analysis only as a means of minimizing intra- 
halo light production by dynamical considerations alone. 
We note that the result of this investigation, despite the 
above caveat, is only a systematically mild reduction in 
IHL across the full spectrum of hos t mass. Our rejectio n 
of Case 2 aligns with the findings of lConrov et alj (|2007l ) , 
in which the authors use a numerically-motivated model 
for the construction of massive galaxies and find that the 
large majority of centrally-merging stellar mass (> 80%) 
must be ejected into the intracluster medium in order to 
reproduce the observed evolution of these central galax- 
ies at low redshift. 

The most obvious trend in Figure 0] is that the IHL 
fraction rises with halo mass from galaxy to group mass 
scales. On average the diffuse fraction is predicted to 
be negligible in M < 10 11 M Q halos and quite substan- 
tial in groups and clusters. This is independent of the 
method used to assign stripped stellar mass to the diffuse 
component or the central object. In both cases, the rela- 
tion between /ihl and halo mass flattens considerably at 
masses above the group scale, tending towards a weaker 
evolution from a diffuse stellar mass fraction of about 
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Gonzalez et al. (2007). The predicted and observed trends show 
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/ihl ~ 20% at a host mass of ~ 10 14 M Q , to a value of 
nearly /ihl ~ 30% at M hos t ~ 1O 15 M . We also see from 
the figure that the initial IHL fraction (without the inclu- 
sion of pre-processed diffuse stellar material) is virtually 
flat on cluster scales, implying that the a priori presence 
of subhalo IHL is largely responsible for the weak increase 
in the host's total diffuse light on those scales. This mild 
trend has also been recovered by numerical simulations 
dMurante et alJl2004t [Monaco et aHl2006tlMurante et afl 
l2007f) as well as observations of intr acluster luminosity a s 
a function of cluster richness (e.g., IZibetti et alj [2005). 

In accordance with our simple model of the previous 
section, the trend of increasing IHL with halo mass is set 
primarily by the convolution of the distribution of subha- 
los that are disrupted (Figure [1]) with the mass-to-light 
ratios of halos (Figure . Consequently, the trends pre- 
dicted by our full model follow closely our general expec- 
tations described in § [3] Specifically, galaxy halos with 
M ~ M t ~ 5 x 10 11 Af Q have massive central galaxies 
because they sit in the valley of the M/L c curve; how- 
ever, these galaxies have low diffuse light contributions 
because they accrete and destroy most of their mass in 
subhalos of mass M ~ 2.5 x 1O 1O M0, where star forma- 
tion is inefficient. Halos at the group scale (~ 10 13 M Q ) 
accrete large numbers of subhalos near the valley of the 
mass-to-light ratio curve. These accreted satellites are 
a copious source of stellar material for diffuse light in 
groups. The diffuse light fraction begins to flatten above 
the group scale because both the host and destroyed sub- 
halos have masses M > Alt, which corresponds to a 
regime where the M/L c relation follows an approximate 
power law. In this case, the ratio of destroyed satellite 
luminosity to central host luminosity is independent of 
mass, L c (M sat ) / L c (Mh ost ) ~ constant. We note here 
that there is also a subdominant effect that contributes 
to the flattening of /ihl at high masses, namely that 
more massive host systems typically accrete their mate- 
rial more recently. This leaves relatively little time to 
disrupt satellites (sec Z05) and results in a lower fraction 
of diffuse, stripped material. 

In order to explore how the total stellar mass within 
halos is divided among the various components (IHL, 
BHGs, satellites) and to more directly compare our re- 
sults with the variety of observational estimates in the 
literature, Figures [5] and [5] show two alternative quan- 
tities. In Figure [5l we ignore surviving satellite galax- 
ies altogether in order to determine the relative impor- 
tance of IHL as compared to the total stellar mass in 
the BHG+IHL. In Case 1 (left) the IHL dominates the 
BHG on cluster scales, contributing 80 — 90% of the com- 
bined stellar mass, while the fraction IHL/ (IHL+BHG) 
declines to ~ 1% on galaxies scales, where it is nearly 
identical to our definition of the intrahalo light frac- 
tion. Again, this is easy to understand in terms of the 
empirically-determined mass-to-light ratios of halos in 
hierarchical cold dark matter cosmologies. On galaxy 
scales, the host halo forms stars at near maximal effi- 
ciency, while its accreted substructures carry compara- 
tively little stellar mass. The IHL/(IHL+BHG) fraction 
is nearly equal to /ihl because nearly all of the luminos- 
ity in non-diffuse (or, for that matter, diffuse) form is 
in the BHG. As host halo mass increases, the efficiency 
of galaxy formation in the central system itself declines, 
meaning relatively more of the non-diffuse light is carried 



10 



0.1 



0.01 



3. 0.1 



t5 



0.01 



0.1 



0.01 



: 1 1 1 1 1 1 1 1 : 

// /■'" ' "'M 

r ' \ 
V / \ 

- / 1 - 
! ' 1 ' 1 ' 1 ' 1 ': 


: 1 1 1 1 1 1 1 1 : 

r v \ 
r m 

// 

- // / - 
': ■ " 1 f ■ 1 ■ 1 ■ 1 ; 


//.' \ 

ft ! 1 
r / M hosl =io' 3 m s ■ 

: ' 1 Li- 1 ' 1 ' 1 : 


//''' \ 

■ /v 

...//..... 


/ ■ 

/ M 

// ; 

// M hosl =10 15 M 
'. ■ i ■ i_Z L ■ i ■ i ! 


10 6 10 8 10 10 10 12 

K ( M ©) 

oil occreted galaxies 

disrupted galaxies 

surviving galaxies 



10 4 



10 e 



10 8 10 10 
K ( M o) 



10 



12 



Fig. 7. — Similar to Figure[T] the differential contribution to the stellar mass from subhalos with an initial stellar mass M* as a function 
of M*. We show this differential contribution for five host halo masses indicated in the legends of the figure. Red dashed lines indicate 
stellar material contributed by disrupted satellite galaxies, while green dot-dashed lines indicate the fraction of stellar mass in surviving 
subhalos for a particular host. 



by the satellites that are not shredded. This causes the 
IHL/(IHL+BHG) fraction to increase more rapidly with 
mass than /ihl- Importantly our result compares favor- 
ably to the ~ 80% IH L to IHL+BHG fraction found by 
iGonzalez et al.l (|2005f) in galaxy clusters. 

Figure [6] depicts a related quantity the IHL + BHG 
fraction relative to the total stellar mass. The IHL + 
BHG fraction is anti-correlated with host mass, decreas- 
ing from ~ 40% on group scales to ~ 30% within large 
clusters. The trend follows from the same logic used in 
the previous paragraph. In addition to the evolution of 
/ihl with mass, the BHG becomes increasingly less lumi- 
nous relative to the sum of the luminosities of its satellite 
galaxies as halo mass increases. The open points with 
error bars show the same quantity deriv ed observation- 
ally fo r individual clusters and groups by Gonzale z et al.l 
(|2007f ) 6 . The predicted and observed trends are remark- 

6 The data table was kindly provided by A. Gonzalez. It included 
information on halo masses, M500, within a radius, -RsoOi where the 
overdensity is 500. The table also included measured IHL + BHG 
fractions within -R500 and within -R200 ~ the radius corresponding 
to an overdensity of 200. For the points plotted on Figure [6] we 



ably consistent, especially on average. Given the observa- 
tional uncertainties, the variance in the observed points 
at fixed mass is also consistent with our prediction, how- 
ever there is a tendency for the data points to skew into 
the upper range of our model's scatter. This may reflect 
a bias in the observational sample, which is selected to 
include systems with dominant BHGs. Indeed, a positive 
trend between dominance of the central BHG and IHL 
fraction is seen in our models (see § 14. 2\ . 

Figure [7] shows the average fraction of diffuse light that 
comes from satellite galaxies of a given stellar mass M„, 
for several choices of host dark matter halo mass. We see 
that the diffuse component (or stellar halo) around small 
Mhost ~ 10 n M Q {e.g., M33) dark matter halos is built 
up from disrupted dwarf spheroidal-type galaxies with 
M* ~ 1O 6 M . Stellar halos around larger Milky- Way- 
converted the reported masses to our A„; r = 337 convention for 
virial mass and plotted the average of the IHL + BHG fractions 
within -R500 and within ifeoo- The error bars reflect the larger of 
the two reported measurements. These corrections amounted to 
~ 20% and ~ 5% changes in mass and (IHL + BHG) fractions, 
respectively, and do not affect the overall trends in any significant 
way. 
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Fig. 8. — The median diffuse light fraction (Case 1) at fixed host mass as a function of the number of massive satellites surviving over the 
halo's history. In the left panel we focus on 10 12 Mq halos (median /ihl ~ 1.7%)and look at the number of surviving subhalos with more 
stellar mass than M* = 10 9 (diamonds) and M* = 10 8 Mq (squares). In the right panel we consider a more massive. 1O 13 - 5 M host halo 
and show the diffuse light fraction (median /ihl ~ 16%) as a function of surviving subhalos with more stellar mass than M, = 10 10 M Q 
(diamonds) and Af* = 10 95 (squares). Though not shown here, the 68% scatter about the median for the 1O 12 M0 host is roughly constant 
at ~ ± 0.3 in logio, while the 10 13 ' 5 Mq host exhibits a smaller variance of ~ ± 0.1 — 0.2 in logio that grows slightly as N su tv increases. 

''- 1 ' " r " '•-• ;h i'-,.i>, .i..-..-r survivors within the groups. Our analysis indicates that 

for galaxy-sized halos, the 68% scatter in each N surv bin 
differs by roughly a factor of two from the bin's median 
value and is approximately constant across the range of 
Nsurv In group-scale hosts, the variance is generally 
smaller (~ ± 0.1 — 0.2 in logio (N surv )) and increases 
slowly as the number of surviving massive galaxies grows. 

This result may explain why some of the 
iGonzalez et al.l (|2007|) clusters have higher IHL fractions 
than we predict (e.g. Figure [6]). These clusters were se- 
lected to have clearly dominant BHGs - in other words, 
to have a less dominant bright satellite population. 
Based on Figure [H we would expect these systems to 
have higher IHL fractions than typical clusters of the 
same mass. 

4.3. Tests for Robustness 

To be sure, our model has several uncertain and poorly 
constrained elements. Particular examples include the 
criterion for removing light from bound satellites and as- 
signing it to the IHL, as well as the evolution of stellar 
mass with time. Our argument in § [3] indicates that 
the overall trends for the IHL that we describe are set 
primarily by the convolution of the mass function of ac- 
creted subhalos with the mass-to-light ratios of accreted 
objects as a function of mass. Further, because we are 
interested in the IHL fraction relative to the total lumi- 
nosity of the system, errors in the normalization of the 
stellar mass function tend to offset each other, if not di- 
vide out precisely. These lines of reasoning suggest that 
the IHL trends that we outline should be robust, at least 
at the qualitative level, but likely at the quantitative 
level as well. Nevertheless, we have subjected our model 
to significant variations in parameter values to assess the 
robustness of the fiducial result. 
Our free parameters govern 

• Dark halo and galaxy disruption via f cr it ~ the frac- 
tion of the initial halo circular velocity that defines 
the critical circular velocity below which satellite 
galaxies are deemed "disrupted" and their stellar 
mass is added to the diffuse light; 



type galaxies, Mh os t ~ 10 M Q , are built from dwarf- 
irregular-size systems, M* ~ 10 8 5 M Q , and intracluster 
light is produced by ma ssive galaxies, M* ~ 1O 11 M 
(see iMurante et al.l 120071 for a similar result from nu- 
merical simulations of intracluster stars). This fact is 
likely to be an important ingredient in understanding the 
metallicities of diffuse stellar components as a function 
of ga laxy luminosity (jMouhcine et al.l 120051 : iFergusor] 
|2007|) ; specifically, more luminous galaxies are expected 
to be surrounded by more metal-rich stellar halos be- 
cause their halos are formed from more massive satel- 
lites. Additionally, note that the differential stellar mass 
distributions become more sharply peaked as host halo 
mass increases from galaxies to groups, reflecting the in- 
crease in relative subhalo luminosity as we approach the 
M ~ M t valley in the M/L c relation (Figure 0). Cor- 
respondingly, the distributions broaden once more as we 
consider the most massive hosts because their subhalo 
populations have moved in large part to the right of the 
valley. 

4.2. The Distribution of the Diffuse Light Fraction at 
Fixed Halo Mass 

A second important feature of the diffuse stellar mass 
fraction is the relative scatter at fixed mass, particularly 
in low-mass halos. The width of the distribution is driven 
primarily by differences in mass accretion histories of ob- 
jects of fixed final halo mass, including the stochastically- 
driven properties of the host's recent merger events and 
the particular orbital parameters for each plunging satel- 
lite. As a general rule, we expect halos that acquired 
their mass more recently to have had relatively less time 
to disrupt the subhalos they host and to have less IHL, 
while early-forming host halos will display the opposite 
behavior. Continuing with this logic, the number of 
bound satellite galaxies should anti-correlate with the 
IHL fraction in objects. Indeed, Figure [5] illustrates that 
our model predicts just such an anti-correlation between 
satellite galaxy abundance and IHL fraction. Of partic- 
ular note is the tight correlation that emerges for group- 
scale objects when considering only the brightest of the 
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• Star formation via a - the star formation parame- 
ter defined in Equation 

• Galaxy luminosities via M/L c as a function of mass 
- the function set from large-scale galactic observa- 
tions to relate central galaxy luminosities and dark 
halo masses. 

In Figure [9] we plot the median IHL fraction computed 
for our fiducial model parameters (thick solid) along with 
various other choices. Note that for simplicity we have 
neglected the "pre-processed" IHL contribution in this 
set of tests. We find that changing the star formation 
parameter a over a very wide range (0.5 < a < 2) pro- 
duces global IHL trends that differ by less than a factor 
of two from the fiducial case. Predictably, the choice of 
fcrit = 0.15 with fiducial star formation (a = 1) results 
in less diffuse light across the full mass range because a 
dark matter subhalo is required to be more severely af- 
fected by the host potential before relinquishing its mass 
to the IHL. However, even this drastic adjustment to f cr it 
produces IHL values that are within a factor of two of 
the fiducial result. 

The most visible change to our main result comes from 
revising our adopted M/L c from the Yang et al. (2003) 
inference to an alternative form advocated by Vale & Os- 
triker (2004) (see Figure 2). The Vale & Ostriker (2004) 
M/L c relation has a steeper "valley" and, as could be ex- 
pected from our discussion in § [31 gives rise to a steeper 
/ihl relation. Even in this case, the overall increase in 
IHL is no greater than a factor of ~ 2 at the cluster scale, 
while the steep faint-end slope of the Vale & Ostriker 
M/L c relation suppresses the diffuse light in small galax- 
ies to below fiducial levels. Despite our limited knowl- 
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Fig. 9. — An assessment of the robustness of our model for the 
IHL fraction (Case 1) under the variation in several inputs. The 
thick, solid line shows our fiducial result for the IHL fraction as 
a function of halo mass. The dotted and dash-dotted lines show 
variations in our prescription for the evolution of stellar mass with 
time [Eq. |5J. Specific parameter values are shown in the legend in 
the lower right corner of the plot. The long-dashed line shows a 
variation in the parameter f cr it that describes the point at which 
the luminosity of a subhalo is assigned to the IHL. The dash-triple- 
dotted line represents the Case 1 result using halo concentration 
parameters exactly twice their fiducial values, to imitate the effect 
of cooling baryons in the central regions of a halo. Finally, the 
thin, solid line shows the result of using o ur fiducial model f or the 
IHL along with th e mass- to-light rat io of I Vale & Ostrikerl (120041) 
rather than that of Yang ct al. (2003). All models show the same 
gross features. Quantitatively the model is very sensitive to the 
mass-to-light ratios of infalling objects in the range within which 
this input can be reliably constrained by independent means. 



edge of star formation, the overall trend appears robust. 
The sensitivity of the IHL fraction to the assumed mass- 
to-light ratios for infalling objects over the range within 
which the mass-to-light ratios can be reliably constrained 
suggests that the uncertainty in this ingredient is a fun- 
damental limitation to the quantitative accuracy of any 
study based on this or similar approaches. In particular, 
the IHL in group-sized systems relies on the precise lo- 
cation of the M ~ M t trough in Figure O Alternatively, 
Figure indicates that it may be possible in the future 
to constrain the M/L c relation between small halos and 
central galaxies by measuring the slope of the IHL frac- 
tion as a function of host halo mass, though more accu- 
rate theoretical methods would need to be employed in 
order to bring this goal to fruition. 

We have managed, throughout this paper thus far, to 
avoid concerning ourselves too much with the baryon 
dynamics that play a significant role in the central re- 
gions of the halos that comprise our dark-matter model. 
However, in the innermost kpc of a dark-matter halo, 
the mass density can be dominated by baryonic mate- 
rial, thereby increasing the concentration parameter typ- 
ically assigned to an NFW density profile and making 
subhalos harder to disrupt. The response of the dark- 
matter profile to the presence of these cooling baryons 
is oft en modeled via adiabat i c contraction (AC) (see, 
e.g., iBlumenthal et alJ 119861 : iGnedin et all l2004f ). In 
cosmological simulations incorporating radiative cool- 
ing as well as sta r formation and supernova feedback, 
iRudd et al.l (|2007j ) find that the effects of baryon contrac- 
tion can be approximated by increasing halo concentra- 
tions at fixed mass by a uniform factor, nearly constant 
over the mass regime of relevance and slightly smaller 
than a factor of two. Since this adjustment applies to 
host halos as well as subhalos, the effect is two-fold: satel- 
lites become relatively more resistant to disruption, while 
their hosts have higher central densities and stronger 
tidal fields as a result of this contraction. In order to 
test the robustness of our model to these competing phe- 
nomena, we double the initial concentration parameters 
of each host halo and subhalo. The resultant IHL frac- 
tion, shown in Figure [9l has a slightly shallower slope in 
the "break" between power-law M/L regimes, although 
the galactic-scale diffuse light is slightly larger that of 
the fiducial Case 1, and the intracluster stellar mass de- 
creases by less than a factor of two. Overall, the effect 
is less important than some of the other uncertainties 
we consider here. We interpret the relative changes from 
our fiducial model to be a consequence of the fact that 
the IHL in larger host halos is governed by the disrup- 
tion of subhalos that are relatively small in comparison 
to the host, ~ 10 n - 5 M Q . These satellites will be more 
resistant to tidal disruption if they are contracted. The 
IHL in small hosts is set by the most massive subhalos 
that merge, and these systems are strongly affected by 
dynamical friction. The dynamical friction force in the 
host halo will be enhanced because of the contraction and 
this enhances massive satellite destruction probability. 

In closing, we reiterate that IHL fractions will be natu- 
rally less susceptible to fiducial normalizations, and that 
our intra-halo light predictions are driven primarily by 
the shape of the M/L c function convolved with the accre- 
tion histories of Figure [T] We conclude that our general 
prediction is robust; explicitly stated, that the IHL frac- 
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tion should rise from a very small value < 1% in low-mass 
galaxies to an appreciable fraction > 20% in cluster-sized 
systems. 

5. DISCUSSION 

The main conclusions of our work may be summarized 
as follows: 

1. The IHL fraction in dark matter halos of mass M is 
expected to increase dramatically from ~ 0.5% to 
~ 20% as we examine systems from the size of small 
spiral galaxies (M ~ 10 n M©) to galaxy groups 
(M ~ 1O 13 M ). The IHL-mass relation becomes 
flatter at a value of - 20% for M > 1O 13 M , in- 
creasing weakly thereafter to ~ 30% for host halos 
of mass M ~ 10 15 M Q . While varying the empirical 
mapping between halo mass and galaxy luminosity 
can produce a slightly higher cluster IHL fraction, 
~ 40%, the overall trends are very robust and are 
governed by the well-known fact that galaxy for- 
mation efficiency varies as a function of mass scale 
while dark matter accretion processes are roughly 
self-similar. Specifically, the subhalos that "build" 
galaxy halos have much lower luminous mass frac- 
tions than the subhalos that build galaxy groups. 

2. The IHL component within galaxy halos is domi- 
nated by the disruption of satellites of stellar mass 
~ 1O 8 5 M while the IHL component in clusters is 
built from more massive stellar systems ~ 10 Mq. 
We expect that more massive galaxies will there- 
fore be surrounded by more metal rich stellar ha- 
los, as has been sugg ested by recent observations 
(jMouhcine et al.ll2005L although lFergusonll2007l dis- 
putes this claim). 

3. The variation in IHL fraction from system to sys- 
tem at a fixed halo mass is driven by variations 
in the accretion history. Systems with fewer sur- 
viving satellites tend to have higher diffuse light 
fractions. The scatter at fixed mass is larger in 
galaxy-sized halos because the light tends to be 
dominated by a small number of massive satellite 
accretion events. As indicated by Figure [SJ the 
number of surviving satellite galaxies in a group is 
expected to negatively correlate with that group's 
IHL fraction, providing an observational expecta- 
tion which future surveys may potentially address. 
This phenomenon may also provide insight regard- 
ing th e comparison of our re sults to observation, in 
which Gonzalez et al. ( 2007) finds a slightly higher 
IHL fraction than our model predicts for group- 
scale hosts, possibly due to a selection effect in 
which their sample systems are typically domi- 
nated by their central galaxies, with relatively few 
bright satellites and thus a systematically larger 
IHL value. 

Current observations place loose constraints on the 
diffuse light fraction on every mass scale. By all in- 
dications, IHL accounts for less than a few percent of 
the tot al stellar mass in large galaxy-sized host ha- 
los (see Sic eel et al.1 [20021: [Guhathaku rta et alJ 2005. for 
discussions concerning the Galactic halo and that of 
M31, respectively), while the diffuse stellar components 



of cluster-sized ho sts are typically about one order of 
magnitude h igher (iMihos et al.ll2005t iZibetti et~aLl l2005t 
iKrick et al.ll2006t iGonzalez et al.ll2007h . A pronounced 
"break" i n the diffuse light belo w the cluster scale is even 
reported (jCiardullo et al.ll2~004r ) . These results are in gen- 
eral agreement with our expectations. 

We predict that the diffuse component around small 
spiral galaxies will contain a very small fraction of the 
primary galaxy's light on average, /jhl ^ !%■ It is inter- 
esting to consider the surface brightness limit that may 
be required to observe such a diffuse component. In Fig- 
ure [10] we investigate a simple example case where we 
have distributed all of the diffuse light predicted for a 
low-luminosity galaxy, L c ~ 4 x 10 9 L Q , into an NFW 
halo density profile that mirrors that of the host halo. 
Solid and dashed lines correspond to the median and 95 
percentile predictions. Here, we have assumed a stel- 
lar mass-to-light ratio of 1 in the R-band. For reference 
we also plot the exponential surface brightness profile 
(|Kim et al.ll2004[ ) for the disk of a system of comparable 
luminosity, the Sculptor group galaxy NGC 300, which 
was shown by Bland-Hawt horn et al.l (|2005f ) to extend 
~ 15 kpc from the galaxy's center without revealing any 
underlying diffuse component. According to our analy- 
sis, a survey reaching ~ 17 kpc from the galaxy's cen- 
ter and achieving 32 magnitudes per square arcsecond 
might be able to detect a stellar halo around NGC 300 
if the diffuse component is comparatively bright, while 
a more average IHL value for the system would require 
an even deeper search. Similar analyses for Milky- Way- 
sized stellar halos indicate that the IHL begins to sep- 
arate itself from a (face-on) disk profile at roughly 29- 
32 magnitudes per square arcs e cond, which is in line 
with the results of llrwin et all (|2005h for M31. This 
provides some idea of the observational depth that will, 
in the future, be required to identify remote stellar ha- 
los around small spiral galaxies. It is worth pointing 
out that some fraction of this light may be in the form 
of recently-destroyed satellites, which should produce 
higher-surface brightness features and will be more eas- 
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Fig. 10. — The R-band surface brightness profile as a function 
of radius for the diffuse stellar component in a small host halo 
(L c ~ 4 X 10 9 Lq), where the IHL is assumed to trace the back- 
ground projected-NFW density profile. Shown are the stellar halo 
profiles for the two values at either extreme of this host's IHL 95% 
distribution, as well as the median value of /ihl — 0.003. For com- 
parison, we plot the surface brightness of the exponential disk of 
a similar system, the Sculptor group member NGC 300, a galaxy 
whose disk extends to at least 15 kpc without the detection of an 
underlying diffuse component (Bland-Hawthorn ct al. 2005). Our 
results demonstrate that this is not unexpected. 
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ily seen (e.g., iBullock &; Johnsto"n]|2005f) . Of course, the 
likelihood of a recent accretion will decrease for lower 
mass galaxies. A more detailed investigation of these is- 
sues is warranted, but the complexity inherent in study- 
ing these issues places such an attempt beyond the scope 
of this paper. 

It is worth noting that while we have focused on ac- 
creted material as the source of diffuse light, several other 
sources hav e been discussed. Th ese include in situ star 
formatio n dGerhard et all l2002f). ejectio n from binary 
systems (|Hollev-Bockelmann et all I2005D. dry merger s 
(in clusters) between ellipticals ([Stanghcl lini et al.ll2006j) , 
and collisionless evaporation (|Muccione fc Ciottil [2004) . 
However, as this work demonstrates, galaxy disruption 
provides a reasonable and seemingly inevitable mecha- 
nism for producing IHL on all scales, although these phe- 
nomena almost certainly have some subdominant role in 
the buildup of a host halo's diffuse light. 

According to our picture, the driving force behind the 
creation of IHL on every mass scale is the stellar mass 
spectrum of intrahalo progenitors (Fig. [7]) . The proper- 
ties of a system's diffuse luminous component can be un- 
derstood as the result of the stochastically-driven merger 
history of stellar-rich satellite galaxies, indicating that 
future observations of intrahalo light could be used as a 
probe of a galaxy's merger history. The predicted trend 
with IHL fraction and halo mass is certainly within the 
scope of future observational work. Interestingly, prelim- 
inary results from the Galaxy Halos, Outer disks, Sub- 
stru cture, Thick disks, a nd Star clusters (GHOSTS) sur- 
vey (|de Jong et al.l -2007) suggest that the stellar halos of 
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low-mass galaxies are, indeed, less prominent than those 
of more massive galaxies. Ongoing surveys of this type 
will be able to test whether the expected trend carries 
over to other mass regimes. As we demonstrated in Fig- 
ure [9j while the qualitative positive trend between IHL 
fraction and halo mass is robust, the slope of the relation 
is sensitive to the underlying relationship between halo 
virial mass and galaxy luminosity on dwarf galaxy scales. 
In principle, this measurement, in concert with models 
of the kind we present, will help constrain the nature of 
galaxy formation in dwarf -irregular- size halos and test 
the accretion histories of dark halos on small scales. 
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